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THE  COUEITE  FLOW  BETVEEN  TWO  PARALLEL  PLATES 
AS  A FUNCTION  OF  THE  KNUDSEN  NUMBER 


I.  INTRODUCTION 

In  a previous  rerxjrt^  a formal  solution  vas  obtained  for  the  heat 
transport  througli  a gas  between  two  parallel  plates  as  a function  of  the  Knudsen 
number  d/\  (d  Is  the  distance  between  the  plates,  X.  is  the  mean  free  path  of  the 
molecules).  The  st£irting  point  was  the  linearized  Boltzmann  eijuation,  since  it 
was  assumed  that  the  magnitude  of  the  disturbance  from  equilibrium,  measured  bjr 
the  retio  AT/T  (lAT  •=  temperature  difference  between  the  plates,  T = average  tem- 
jKratiire),  was  small. 

In  the  following,  the  same  method  will  be  applied  to  the  problem  of  the 
Couette  flow  of  a gas  between  two  parallel  plates  as  a function  of  the  Knudsen 
number  d/X.  We  will  again  assume  that  the  Mach  number  which  in  this  case  is  the 
ratio  of  the  average  flow  velocity  to  the  mean  molecular  velocity,  is  small,  so 
that  the  disturbance  of  the  equilibrium  due  to  the  moving  plate  is  also  small  in 
this  case. 

Since  the  method  and  the  general  features  of  the  solution  are  quite 
similar  to  those  of  the  heat -transport  problem,  only  an  outline  of  the  calcula- 
tions will  be  presented  in  sections  IT,  III,  and  P/.  In  section  V some  diver- 
gence difficulties  will  be  discussed  which  also  occur  in  the  heat -transport 
problem,  and  which  are  due  to  the  parallel-plate  geotautry. 


. S.  Wang  Chang  and  G.  E.  Uhlenbeck,  "’file  neat  Transport  between  Two  ParaU.el 
Plates  as  Functions  of  the  Knudsen  I'dimber”,  Univ.  of  Mich.,  Eng  Res.  Inst., 
Project  M999,  Sept.,  19^3. 
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II.  FOEMjLATION  OF  THE  PROBLEM 


We  t.ake  the  y-z  plajie  halfvay  "between  the  plates.  The  upper  plate, 
situated  at  x - d/2,  is  stationary,  while  the  lower  plate  at  x = d/2  is  moving 
In  the  z direction  with  a velocity  w (measured  in  units  of  *v/m/2kT^) . The  velocity, 
w,  is  sissumed  to  he  much  smaller  than  unity  so  that  terms  of  the  order  w2  and 
higher  will  be  neglected.  The  notation  of  the  previous  report  will  he  used. 

The  distribution  function  is  written  as 

J 1 + ^ ^ . vvr)] 

It  is  convenient  to  take  for  the  zeroth  approximation  distribution  function 


where  C ■ c - w/O,  n is  the  equilibrium-number  density,  and  T the  equilibrium 
temperature,  which  are  all  constants.  The  Boltzmann  equation  is  again 

where  J is  the  collisior  operator. 

To  formulate  the  bo\indary  conditions  we  introduce  the  distribution 
functions  for  the  molecules  going  up  and  down,  i.e.,  f'*’  and  f",  and  the  corre- 
sponding disturbances  hC  and  h”.  In  terms  of  the  h's,  the  bounlary  conditions 
are 

V(~  z)  = ^ + ( '"*^0  ("  4 . " Cx)  X > O (io.3 


r:  o<  ( g)  - wCp  -V  C l-0()  4 ) ~ X < o Cit} 

where  Ct  is  the  accommodation  coefficient  and  the  constants  B'*’  and  B~  are  to  be 
determined  by  the  conditions 
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expressing  the  conditions  that  the  number  of  molecules  per  unit  area  between  the 
plates  is  given  (equal  to  nd)  end  that  in  the  steady  state  there  is  no  net  flow 
In  the  x-direction.  Up  to  terms  linear  in  w we  have  the  symmetry  condition 

and  consequently 

y ) ^ - C-Cx,  -xj)  , 

With  the  symmetry  condition  (4),  Eq.  (3a)  is  automatically  satisfied,  and  the 
boundary  conditions  (2a)  and  (2b)  oire  equivalent  with 

= - e>“  - 6 . 

Thus  our  problem  reduces  to  the  solution  of  the  integral  differential 
equation  (1),  subjected  to  the  boundary  condition  (2a)  and  the  axuclllary  condi- 
tions (3b)  and  (4).  The  physicaJ.  quantities  we  are  interested  in  are  the  drag 
on  the  upper  plate,  the  velocity  distribution  as  a function  of  x,  and  the 
velocity  Jumps  at  the  plates. 


III.  GENERAL  SOLUTION 


At  first  glance,  this  problem  seems  to  be  more  complicated  than  the 
problem  of  the  heat  conduction  b*?cause  of  the  lack  of  axial  symmetry  around  the 
x-axis  in  the  present  problem.  It  will  be  seen,  however,  that  this  is  not  a real 
dlfflcxilty.  It  is  true,  though,  that  in  the  development  in  eigenfunctions,  we 
must  now  use  the  spherical  hai’monlcs  instead  of  the  Legendre  polynomials  used  In 
Ref.  (1).  We  write 
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where  8 is  the  angle  between  the  velocity  C and  the  x-axls  and  0 1b  the  azimuthal 
angle  meesured  from  the  y-axls,  and  develop  h according  to 

'I  L 

* — 

Because  of  the  reality  of  h one  must  have  aj./m  " ®Tim’  I'Tom  the  symmetry  of  the 
pi’oblem  with  regard  to  the  y direction,  it  follows  that  h must  be  an  even  fiinction 
of  Cy,  which  has  as  a consequence  that  the  a^^j^  for  odd  m must  be  pure  imaginary 
The  first  few  of  the  development  coefficients  Sixlm.  related  to  quantities  with 
physical  interest.  For  Instance,  the  density  and  temjjerature  are  found  from 


n *'4  - 7 

T(«)-  T[t-f 


and  the  physical  quantities  (shearing  stress,  average  velocity,  and  heat  flux, 
all  in  the  z-dlrection)  in  which  we  are  especially  interested  are  given  by 


3 o 


^ ®-9l  1 

^ KToii 


(W 

(feb) 

v.Gc  ^ 


From  the  conservation  laws  it  follows  tha*:  c^,,  Pxx>  Pxy>  Pxz>  'lx  '^on- 
stantc . We  will  require  that  Cj,  - 0 (see  eq.  (5b));  p is  zero  by  symmetry. 
This  implies  that 


^ 01  c ~ 

0 

^ 0 J.I  ^ 

C/huiA". 
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^«oo 
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*There  are  some  misprints  in  the  corresponding  formulas  in  the  previous  report, 
Ref.  (1) 


k 


ENGINEERING  RESEARCH  INSTiTUTE 


UNIVERSiri  OF  MICHIGAN 


¥ 


t- 

k 


i 

f 


s. 

¥• 

'r 

t 


L. 


Introducing  the  development  of  h into  the  Boltzmann  equation,  we  find 
that  our  prohlem  reduces  to  the  solution  of  an  infinite  set  of  homogeneous  linear 
differential  equations : 


=.  -n  J a [ _L  vt  ,r/  '1 


with  the  "boundary  conditions : 


q) 


. e\-<> 


Finally,  = 0,  and  from  the  synmetry  property  eq.  (4)  it  follows  that 


.even 


.odd 


^Im 


1 


in  X according  as  1 - m isj 


odd 


even. 


•The  square  bracket  in  eq.  (7)  is  different  from  zero  only  when  m = m',  and 
1 - i'  is  odd.  Thus  coefficients  with  different  volxies  of  m are  not  coupled. 


Since  forlral>l  the  integral  in  the  boundary  condition  (8)  is  zero,  it 
is  clear  that  the  equations  for  with  lm|>l  are  compl->itely  homogeneous,  so 

that  the  only  solution  will  be  zero.  This  is  nlso  the  case  for  m = 0,  because 
in  eq.  (8)  the  qixantlty  B is  also  an  unknown  Euid  the  inhomogeneous  peurt  is  only 
the  term  containing  wC^.  This  term  is  also  zero  for  m = 0.  Hence  one  may 
conclude* 


b = o 


a 


y1-«vi  = o , for  m / + 1 


We  use  the  same  procedure  eis  was  employed  in  Ref.  (l).  One  must 
distinguish  between  even  and  odd  values  of  £.  Eliminating  \f2t  + 1)1 
tains  an  infinite  set  of  second-order  homogeneous  differential  equations  for 
®r(2l)l»  which,  with  the  use  of  the  symmetry  conditions,  has  as  solutions: 


I 0 
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L._ 


*This  Is  In  contrast  to  the  case  of  heat  conduction.  It  is  also  physically 
evident  that  a temperature  .rpradlent  will  produce  a change  of  the  density  n near 
the  pla.s,  while  this  is  not  the  case  for  a velocity  gradient. 
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where  the  P^'s  are  the  positive  and  nonzero  solutions  of  the  secular  determinant 


A = 


= o 


and  where  in  the  sets  of  the  ^j.(2i)l  ^ taken  as  arbitrary,  for  which 
we  choose  hof\  and  . Furthermore,  from  the  constaincy  of  it  follows  that 
hoai  ■ 0 for  aH  1. 

The  coefficients  + i)i  completely  t’etermined  and  one  finds 


\ --  I c'pi t»- »,  'Vv'itt's  ] ^'r'2.K'\ 

The  remaining  unknown  constants  ho°i  and  have  to  be  determined 

from  the  boundary  conditions.  For  this,  one  needs  first  the  a^(2i)l»  which 
one  obtains; 

where 

Putting  eq.  (11)  in  the  boundary  conditions,  one  ib  led  to  the  set  of  inhomogene- 
ous linear  equations 


Cli) 


where  the  known  constant  matrices  B are  defined  by 

^ u'' 
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In  the  same  sense  as  iu  Kef.  il),  eq.  (12)  is  a set  of  Inhomogeneous  linear 
equations  with  the  same  number  of  equations  as  that  of  unknowns.  Thus  the  prob- 
lem is  formally  solved. 

The  complete  solution  involves  infinite  determinants.  We  are,  as  yet, 
unable  to  discuss  the  convergence  of  these  determinants  nor  the  convergence  of 
the  "breaking  off"  processes  which  we  will  use  &b  in  Ref.  (1). 


K.  ^ MAXWELL  MOLECULES 

We  will  use  the  same  "successive  approximation"  method  as  in  Ref.  (1). 
Furthermore,  ve  will  restrict  ourselves  to  the  Jfaocwell  molecules  so  that 


0)  . 

and 

•w- 

(^3)  J ~ ® unless  both  r and  I a:^  zero  which  has 

as  a consequence  that  all  zero  except  bdai. 

Since  we  do  not  expect  to  get  anything  from  the  zeroth  approximation,  we  start 
with  the  first  approximation. 

A.  First  Approxiination 

A.S  in  the  heat-conduction  and  sound-propagation^  problems,  we  will  use 
in  the  first  approximation  the  eigenfunctions  for  which  2r  + i 4 5,  and  l<5. 
since  only  the  '|('s  with  m = 1 enter  the  pi’oblem,  we  therefore  need  to  0;rly 

■toil,  and  +021-  We  know  that 

(o) 

^21  - const.  = ho2i-  (i3ct^ 


®C.  S.  Wang  Ctiang  and  G,  E.  Uhlenbeck,  "On  the  Propagation  of  Sound  in  Monatomic 
Gases",  Unlv.  of  Mich.,  Eng.  Res.  IriSt.,  ProJ . K999^  Oct.,  1952. 
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The  other  tvo  a's  sxe 


'>'Wm  [^'vTu  ] , 

~ 'bk>ovv  L"^  "luv  ) 1 . 

C\2c.) 

Since  p = 0 in  this  case,  the  houndary  condition  (12)  reduces  to 


V- 

Oi\  " 


A-o<  , _ "nA  -r-'^  V 

■a/y  ”x  '«>■»>,  e».\ 


Where  j 

Using  the  values  of  S and  T*as  given  in  tables  in  Appendix  II,  one  finds 

V,^°^  = - -T ^ — 

oi-''  A i-x  g\^  X 

‘ a ~ 

where  ' • Using  e^.  (6)  and  (13),  one  obtains  the  following 


C ~ '^V'T  ^ ^ , 

3 X ^ 

7T  _ f i_  1 ] 

h.~  T at  J 

i-°<  3 dL  -> 


('-Vo,) 


(\4t> 


H T'X 


C'4-c) 


r {r- 
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where  w is  still  In  vuiits  s^T/m.  Therefore  the  velocity  slip  at  the  upper 
plate  X = d/2  is 


n 


V - i ~ 

X. 


I = yi-  hzjL.  ^L_ 


i 


2>  ^ A 

^ 3 ir 


For  P)i»^T ^ x/d  and  a = 1,  this  is  in  agreement  with  the  result  of  Maxwell. 


B.  Second  Approximation  (2r  + ^ and 

The  eigenfunctions  entering  the  problem  are  toii^  iiiij  +031?  and 

+151.  The  +'s  with  odd  values  of  I - m are  tosi  ana  Yisi.  The  only  nonzero 
value  of  p is 


The  constants  to  be  determined  by  the  boundary  conditions  are  bo*^!  and  bf§\. 
Equations  (12)  become  a set  of  two  equations; 


and 


C - b,., 


Solving  these  equations,  one  then  can  calculate  the  quantities  of  physical  in- 
terest. It  turns  out  that  in  all  approximations  the  drag,  the  velocity  distribu- 
tion, and  the  velocity  Jump  have  the  form 


"i" -a  - 


^ 4 A_  . 

\ M 3 B 


r . 4r  2.  "E^C  1 

r - Vjl  1 1 - ii.  ii 1-^  A_k I I / > r) 

" [ ^ \ d 2.-<X  XNlT  A r-o(  Xfir  A J . ' ^ 


1 ' ^ 


£ 2^  XrniA-  ^ 4 i-o(  S'Tir  A 

'A  ^.«^s^5  T 2.o(siS 


■2.  3 
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where  the  svm  is  over  all  the  roots  pi,  and  = sinh  (npj^d/2).  The  functions 
A,  B,  and  are  all  functions  of  K = d/\. 

In  this  approxteation,  one  finds 

k - -V  t 

^ io<  ^ 4.0W  ^ > 

icA  tt  ^ > 

c _ BtiiL  -L 

fc-  — \<r  r 

• W ^ , 

where  t •=  tanh  (npd/2),  c = cosh  (npd/2). 

C . Third  Approximation  (2r  + £ 5»  i.  ^ 5) 

The  eigenfunctions  to  he  taken  in  this  case  are  toii»  iosi#  tin*  Yo3i» 
♦211*  +131*  There  are  four  \|f's  with  odd  values  of  f - m.  Hence 

there  are  two  pairs  of  p differing  from  zero.  They  are 

Pi  = 0.7151  Ag  and  ps  * I.5U2  A^  , 

where  A?  = A*  -/x/kT'.  We  find,  further. 


I 


(1) 


(z) 


ho7i  = -0.2^56  h , ht^i  . 2.559  hi'ri. 

Solving  the  set  of  three  equations  for  the  boundary  conditions,  one  finds  then, 

+ ^^(o.t3\8T,4o.s-3?‘Vt»)-v  C.35‘^zt,t2.  ^ 

^ 0.1-]  14  ^ 

P = c. 3000  + 0. nit t.) 

V.  2^  ' 
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D.  Fourth  Approximation  (2r  + ^ 6 , and  £ < §) 

We  T3»*oceed  aa  before,  with  the  four  values  of  p's  and  the  corresponding 
ratios  of  given  In  the  following  table. 


O.U60HA4 


1.065Ai 


I (.  0.711)  5.525) 


l,26lAi 


1.6lOAi 


(3.322)  — (15.58) 


^ (-  1.218) 


(1.362) 


^ (-  6.204) 


^ (1.213) 


J2  (8.442) 


1 ^ j 1 fS 

4^il  ("5.349)  I }7j  Tj[  (32.51) 


'Daklng  a * 1 to  siiqjlify  the  congjutation  we  aoTived  at  the  following  results; 

= a. 3^  + I C.3  V,+  J t^+  l.SO  -t  2.70t,t^t 

i a.73t,t<^  + i-n -tjtj  1 2.8-1  t 

-V  3.03  !•  S'©  + 2. 74  + 3.\0 

B = 1.5?'+  1-4^,+  l-33tv+  l.4t+,-l:^+  l-SS-+,-^  + 

+ 2-4Z-t,-V4t  i.‘)3+i+<v+  + 

i.4ott+^  + i.+r^r,-t-3n  + 2.44,  ^ 

E ( c.s-6:l  + 0-4S7  ti.+  0 S-Pit  t o.4?'?t.^.  + 0.4S3t^t^+ d.4l7t,uf 

-V  o.4?2+^t^  + 0 433  ^ 

(32x]  + 3.i'4t,t  3.3x3rt3+  3,o3'i+<|.+  3.oll+,t^  + 

3.\y3t,-tit^+  3.\iS-4,t3t4)  ^ 
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^3  ~ ■ ( q-"o.S~  -V  o.gi'y  o.  o.  ?53-t-4.T 


-V  0.700  + fl.w  ^ 


l.^^V  ^•7oStT4+  i <V'^S't,t^+ 


^ 3.  + a.sg:'^ 


tor  the  discussion  of  these  results  we  have  plotted  three  sets  of 
cxurves.  Figure  1 contains  a set  of  four  curves  for  Khudsen  afiai^s* 

the  Kiudsen  numiber,  d/\.  Curve  I is  for  the  Stokes-Havier  approximation.  Curves 
II.  Ill,  and  IV  are  results  from  our  first,  second,  and  third  approximations 
respectivelj'.  We  have  not  plotted  the  curve  for  the  foiirth  approximation  "because 
we  do  not  expect  anything  new  and  the  numerical  computation  is  very  laborious. 

The  straight  line  on  the  left  is  the  initial  slope  for  the  exact  solution.  All 
the  curves  from  the  different  approximations  have  the  same  value  (unity)  for 
K = 0,  and  for  K equal  to  infinity  they  all  approach  the  same  limit,  zero.  The 
initial  elope  for  the  four  approximations  are  listed  below. 


Exact 

First  approximation 
Second  " 

Third  *' 


InitlecL  slope 
-1.2U2 
-0.1;25 
-O.I458 

-O.55U 

-0.598 


Figure  2 consists  of  a set  of  three  curves  for  the  first  three  approxi- j 
mations  for  the  velocity  slip  at  the  upper  plate  as  a function  of  K.  The  initial! 
slo’Te  is  given  by;  j 


s 

t#-  • 


First  approximation 
Second  " 

Third 

Fourth 


Initial  slope 

-o.U9it 

-0.6UU 

-0.755 


We  have  drawn  Fig.  5 to  show  the  velocity  distributions  as  functioiis 
of  X for  K = 10.  Since  the  velocity  distributions  for  the  different  approxi- 
mations deviate  very  slightly  from  the  stralglit-line  distribution  of  the  first 
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approxliDifiltJ.on  ve  liave  drawn  instead  of  c^, 
R = 1 - 2Cj,(x)/w  and  R^^  is  the  value  of  R 


against  x,  tlie  curves  Rj^  - R,  where 
from  the  first  approximation,  namely. 


K_  : 

~ 4 ^ X 

•3  4 

As  in  the  case  of  heat  conduction,  the  boundary  effect  is  shown  hy  the  sharp  rise 
of  the  cTjrve  near  the  wall. 


V,  THE  APmOAGH  TO  TEE  KNUPSEN  LIMIT;  A DIVERGENCE  DIFFICULTY 


In  Ref.  (1)  we  pointed  out  that  the  behavior  of  the  exact  solution  for 
the  heat -conduction  problem  is  for  large  K quite  di+'ferent  from  the  behavior  at 
small  Khudsen  niimber.  T^e  approach  to  the  Ciaucivc  gas  limit  (K»l)  is  compli- 
cated by  the  occurrence  of  the  hj’perbolic  functions,  so  that  a development  in 
inv»irse  pxjwers  of  K is  not  possible . He  see  from  section  III  that  the  same  Is  true 
for  the  Couette  flow.  On  the  other  hand,  it  seems  that  the  approach  to  the 
lOmdsen  gas  regime  (K«L)  is  quite  regular  so  that  a development  o^  all  relevant 
physical  quEuititieB  in  powers  of  K should  be  possible. 

In  fact,  in  Ref.  (1)  we  gave  the  first  two  terms  in  such  a development 
for  the  heat  flux.  ArAlogous  results  can  be  fovind  for  the  drag*.  In  the  zeroth 
approximation  one  finds 
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with  w still  in  units  ofNl2kT/m.  In  the  first  approximation 
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For  a = 1 this  reduces  to  the  expression  given  in  Ref.  (5).®  There  the  square 
bracket  which  is  always  nsgatlve  has  been  evaluated  for  elastic  spheres  and 
MsLXwell  molecules  . The  resiilts  are  as  f ol  lows , 


♦These  results  can  be  found  either  by  the  method  descriled  in  Chap.  Ill  of 
Ref = (l'  or  from  the  general  solution. 

^C.  S.  Wang  and  G.  E.  Uhlenbeck,  "Transport  Fnenoaiena.  in  Very  Dilute  Gases", 
CM  579,  UMH-5-F,  Univ.  of  Mich.,  Eng.  Res.  Inst.,  ProJ . m6oU-6,  Nov.  i>,  19-9- 
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For  elastic  spheres  of  diauie!«r  a, 

Jfl  i4'(t) 

V,  lo)  4-?  <> 

For  Maxwell  molecules, 

v>  1 
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In  this  case  the  resixl'c  obtained  frea  tic  general  solution  Is  expressed  in  the 
eigenvalues  of  the  collisioa  operator  in  the  form 
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The  suasffiatlon  can  he  carried  out  (for  details  see  Appendix  I)  and  confirms  eq,. 
(l6b).  The  integral  vas  evaluated  nTjraerl cally  in  Ref.  (3).  with  the  result 


V') 

V _ 

where  - 'rA^ix/iT  • a = 1,  the  exact  value  of  the  InltleLl  slope  in 

Fig.  1 is  -1.2U2. 

A.  difficulty  appears  when  one  wants  to  calculate  the  second  approxi- 
mation of  the  drag,  or  when  one  wants  to  find  analogous  expansions  for  the 
average  velocity  distrihution  or  for  the  velocity  jump  at  the  plates.  Formally, 
one  finds 

and  in  the  first  approximation, 

U) 
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The  bracket  expressions  are  both  divergent  integrals  because  of  tiie  factor  1/c^. 
The  same  difficulty  occ\irs  in  the  hest-ccaducticn  problem,  although  we  did  not 
notice  it  at  that  time.  The  second  approximation  zo  the  heat  flux  given  in 
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Ref.  (1),  p.  8,  eq.  (19)  is  also  divergent,  and  the  temperature  distribution  in 
the  first  approTlination  for  which  one  can  derive  the  formal  expression 
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is  divergent  Just  as  c 
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The  divergence  of  these  integrals  is  of  a logarithmic  nature  and 
is  always  due  to  the  factor  l/c^.  The  origin  of  these  divergence  difficulties 
is  therefore  clearly  the  parallel-plate  geometry.  MoleculoC  which  are  emitted 
nearly  parallel  to  a plate  will  have  to  transverse  a very  long  path  before 
reachiiig  the  ether  plate,  so  that  the  Khudsen  approximation  will  not  be  valid 
for  these  molecules-  Or  one  can  say,  the  average  Knudsen  n>amuer  for  mcleculeR 
e?nitted  in  all  directiono  Ic  logarithmically  infinite  even  if  d/\«L.  We  have 
found  that  by  taking  two  concentric  spheres  (radii  a s.r.5  b)  instead  of  two 
parallel  plates  all  divergences  disappear.  For  d = a - b«$./2)  (a  + b)  = R, 
quantities  ouch  as  the  temperature  distribution  will  in  first  approximation  in 
K = d/x  contain  terms  proportional  to  3.n  (R/’d),  which  blows  up  for  the  parallel 
plate  geometiry. 

One  may  conclude,  therefore,  that  these  divergence  difficialties  will 
not  affect  any  real  physical  situation,  and  that  the  behavior  of  the  solutions 
near  the  Khudsen  limit  can  still  be  considered  to  be  regiaar,  so  that  a develop- 
ment in  powers  of  K is  possible. 
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APPENDIX  I 


DERIVATION  OT  EQUATION  (l8b)  j’ROM  ECgJATION  (^) 


In  terms  of  the  eigenvalues  of  the  collision  operator,  the  first 

Order  correction  to  the  pi^essure  p in  the  iuiudsen  limit  is  given  hy 

xz 
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The  expression  for  Xraf  is 


1 

+ - (l-Sr»8w)| 

In  the  hraces,  we  see  that  the  second,  third,  and  fourth  terms  sxe  obteinshle 
from  the  first  term  hy  replacing  0 hy  « - 9 , 0,  and  jt/S.  Thus,  substituting 
eq.  (2it)  into  eq.  (2?),  we  see  that  we  need  only  to  calculate  the  sum 

Putting  r -f  I = i,  we  have 

The  leist  sxaa  can  he  done  with  the  help  of  the  relation 
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and  thuB  to  eq.  (l8b). 
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